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PREFACE 


An  attempted  asymptotic  series  solution  in  terms  of  inverse  powers 
of  the  axial  distance  for  the  incompressible  laminar  wake  far  behind 
a  very  slender  cylinder  breaks  down  at  an  early  stage.  In  this  Memo¬ 
randum,  the  PI1C  method  is  applied  to  the  problem  and  yields  a  uniform¬ 
ly  valid  asymptotic  solution.  The  method  of  approach  used  should  be 
of  interest  to  researchers  in  fluid  mechanics  concerned  with  the  solu¬ 
tion  of  singular  perturbation  problems,  while  the  results  should  be 
of  interest  to  those  involved  in  wake  studies. 


SUMMARY 


The  Oseen  solution  for  the  wake  far  behind  a  very  slender  cylin¬ 
der  indicates  that  an  asymptotic  series  solution  in  inverse  powers  of 
the  axial  distance  might  be  possible.  However,  if  an  attempt  is  made 
to  carry  this  out,  the  solution  breaks  down  at  an  early  stage.  This 

occurs  when  the  coefficient  of  the  inverse  second  power  of  the  axial 

_2 

distance  (x  )  exhibits  the  wrong  behavior  at  the  edge  of  the  wake. 

This  difficulty  is  overcome  by  applying  Lighthill's  technique  (the 

P11C  method)  for  rendering  approximate  solutions  uniformly  valid.  The 

2 

correct  series  contains  the  term  0(Xn  x/x  );  it  was  the  omission  of 
this  term  which  caused  the  difficulty  encountered  in  the  original  se¬ 
ries  expansion.  In  addition,  it  is  shown  that  the  next  term  in  the 
_2 

series,  0(x  ),  contains  an  indeterminate  factor.  This  indeterminacy 

is  due  to  the  occurrence  of  eigenfunction  solutions,  but  the  reason 
for  their  existence  is  not  explained.  Earlier  work  on  this  problem, 
however,  indicates  that  this  is  connected  with  the  neglect  of  the  ini¬ 
tial  velocity  profile  in  obtaining  the  asymptotic  solution.  The  solu¬ 
tion  obtained  has  some  relevance  to  the  wake  behind  any  finite  axisym- 
metric  body. 
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I .  INTRODUCTION 


The  asymptotic  solution  for  the  incompressible  axisymmetric  wake 
far  behind  a  long  thin  cylinder  (or  "needle")  was  given  in  Ref.  1. 

(This  reference  also  contains  a  discussion  of  the  applicability  of  the 
"needle"  solution  to  more  general  axisymmetric  bodies.)  The  solution 
was  sought  in  the  form  of  a  series  in  inverse  powers  of  the  axial  dis¬ 
tance.  However,  the  solution  so  obtained  was  not  a  uniformly  valid  one. 
This  first  became  evident  when  the  second  approximation  (corresponding 
to  the  negative  second  power  of  the  axial  distance)  did  not  approach 
zero  exponentially  at  the  edge  of  the  wake.  A  detailed  analysis  of  the 
problem  revealed  that  the  difficulty  lay  in  the  assumed  form  of  the  ex¬ 
pansion  and  indicated  that  logarithmic  terms  had  to  be  included.  In 

2 

particular,  it  was  shown  that  the  term  following  0(l/x)  was  0(£n  x/x  ) 

2 

and  not  0(l/x  )  as  was  initially  assumed  (x  is  axial  distance).  The 

2 

series  solution  was  carried  out  to  the  term  0(4n  x/x  ),  and  it  was  shown 

2 

that  within  the  context  of  the  asymptotic  approach,  the  term  0(l/x  ) 
could  not  be  completely  determined.  This  indeterminacy  was  explained 
as  being  due  to  the  neglect  of  one  boundary  condition  in  obtaining  the 
asymptotic  solution.  This  last  boundary  condition  is  the  initial  veloc¬ 
ity  profile  at  the  base  of  the  cylinder.  The  manifestation  of  the  lost 
condition  in  determining  the  terms  in  the  series  was  in  the  form  of 
eigensolutions ,  that  is,  complementary  solutions  of  the  equations  sat¬ 
isfying  the  boundary  conditions  on  the  axis  and  edge  of  the  wake.  Ar¬ 
bitrary  multiples  of  these  eigensolutions  could  be  added  to  certain  terms 
in  the  expansion,  the  formulation  of  the  problem  being  such  that  there 
was  no  way  of  determining  their  contribution. 

In  the  present  Memorandum,  we  wish  to  reconsider  this  far-wake  prob¬ 
lem  and  demonstrate  how  the  correct  asymptotic  solution  can  be  obtained 

by  using  Lighthill's  technique  for  rendering  approximations  uniformly 

(2) 

valid  (the  PLK  method,  or  "method  of  strained  coordinates"'  In  ad- 

2 

dition  to  correctly  predicting  the  0(4n  x/x  )  term,  the  technique  also 

determines  the  explicit  form  of  the  eigensolution  responsible  for  the 

2 

indeterminacy  in  the  0(l/x  )  term. 
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Applicatlon  of  Lighthill's  technique  to  this  problem  does  yield 

the  solution  more  readily  than  the  analysis  given  in  Ref.  1;  however, 

its  efficiency  is  diminished  by  a  loss  of  insight  into  the  basic  cause 

2 

of  the  indeterminacy  appearing  in  0(l/x  )  and  higher-order  terms. 
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II.  BASIC  EQUATIONS  AND  ASYMPTOTIC  SOLUTION 


The  equations  of  motion  governing  the  development  of  the  axisym- 
metric  incompressible  wake  without  pressure  gradient  can  be  written  as 


d  d 

^  <ru>  +  37  <rv> 


du  ,  du 

U77+  V77 


Id  /  du\ 

7  7F  [T  37 ) 


(Although  the  assumption  of  constant  pressure  would  seem  to  limit  the 
resulting  solution  to  very  slender  bodies,  the  solution  does  have  some 
application  to  the  wake  far  behind  any  axisymmetric  body.  See  Ref.  1 
for  a  detailed  discussion  of  this  point.) 

Here,  x,  r,  u,  and  v  are  nondimens ional  variables  defined  by 


/ue\*  li 

U  )  e  u  *  u 


v  = 


where  x^  and  r^  are  cylindrical  coordinates  (x^  is  the  axial  distance, 
and  r^  the  radial  distance  measured  from  the  axis  of  symmetry),  u^  and 
v^  are  the  corresponding  velocity  components,  U  is  the  free-stream  ve¬ 
locity,  v  is  the  kinematic  viscosity,  and  9  is  the  momentum  thickness. 

Equation  (1)  can  be  satisfied  by  introducing  a  stream  function  i|r , 
defined  by 


1  d* 

i  7r 


1  dt 

r  7x 


Equations  (1)  and  (2)  are  to  be  solved  subject  to  the  boundary 


conditions 


u  -  1 


The  solution  is  most  conveniently  expressed  in  terns  of  a  new  in¬ 
dependent  variable  Q,  defined  by 
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2 

r 

4x 


(6) 


which  replaces  r;  we  also  write  the  stream  function  as 

il  -  xf(x,£)  (7) 

In  Ref.  1  the  initial  approach  to  an  asymptotic  solution  is  made 
by  expanding  f(x,£)  in  the  following  series: 


f  (x,0 


00 


I 

n-0 


fn<» 

n 


(8) 


However,  this  method  fails  when  it  is  found  that  £2(0  approaches 
zero  as  C  “*  ®  algebraically  and  not  exponentially,  as  should  be  ex¬ 
pected  on  physical  grounds.  An  alternate  formulation  was  used  in 
Ref.  1  to  indicate  the  form  of  the  correct  expansion. 

In  this  Memorandum,  following  the  methodology  of  the  FIX  tech 
nique,  we  expand  the  stream  function  and  the  x  coordinate  in  the 
following  series: 


*  -  pf(P,C) 


(9) 


x  ■  P  +  x^P)  +  ...  (10) 


2 

where  Q  ■  r  /4P. 

Substitution  of  Eq.  (9)  into  Eq.  (3)  leads  to  the  following  ex¬ 
pressions  for  u  and  v: 


u 


(ID 
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7  t£  -  «£c  +  V  # 


From  Eq.  (10)  we  find  that 


dx  dx/dP  1  +  x^(P)  +  ... 


1  -  x^P)  +  ... 


assuming  that  x*(P)  «  1  when  x  is  large.  Then 

v  -  -  -  [f  -  Cfc  +  P£p](l  -  *{(&)  +  -•)  (W 

Substituting  Eqs.  (11)  and  (14)  into  the  momentum  equation,  Eq. 
(2),  and  using  Eq.  (13),  we  obtain 

Cfccc +  £cc +  [2  ££cc +  I  eVcc  "2  P£bc£:](1  '  xi(e)  +  ,”>  "  0 


In  terms  of  the  functions  f  (C)»  the  boundary  conditions  are 

n 


£n(0)  -  0 


JZ  £'(c>  -  0 


C  -  0 


£0'(C)  -  2 


f '(C)  -  0 

n 


as  £ 


n^l 


If  we  now  substitute  the  series  expansion  of  f(p,£)  given  in 
Eq.  (9)  into  Eq.  (15)  and  begin  equating  coefficients  of  different 
powers  of  P"1  to  zero,  we  find  that  the  terms  independent  of  P  yield 
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£f  "'  +  t"  +  f  f"  •  0  (17) 

^  o  o  2  o  o 

A  solution  of  this  equation  satisfying  the  boundary  conditions  is 
f  (C)  =  2C,  and  this  term  is  simply  the  uniform  free  stream. 

Before  any  further  equations  can  be  obtained,  we  must  choose  x|(g). 

In  Ref.  1  it  was  found  that  f ^ (C>  8oes  t0  zero  as  C"2  in  the  llmit  of 
£  -»  oOj  and  therefore  the  expansion  was  considered  to  have  broken  down 
at  this  stage.  Through  an  appropriate  selection  of  the  value  of  x^P), 
we  shall  try  to  eliminate  those  terms  in  f 2( 0  which  do  not  decay  ex¬ 
ponentially  near  infinity.  Since  \|/  is  expanded  in  a  series  in  inte¬ 
gral  powers  of  0  \  the  simplest  choice  for  x^(P)  is 

x'(P)  (18> 


where  A  is  a  constant. 

If  the  coefficient  of  g”1  is  now  equated  to  zero,  we  obtain  (as 
suming  fQ(£)  ■  2£) 


£fj"  +  (£  +  l)f'  +  t[  -  0  (19) 

_2 

while  the  coefficient  of  0  yields  the  equation 

d“'+  (C  +  1)£'  +  2£'  +  \  t'2  -  §  (2(f”  +  2£{)  -  0  (20) 

Equation  (19)  is  the  same  as  the  equation  for  f^G)  found  in  Ref. 
1,  and  a  solution  satisfying  the  boundary  conditions  is 

£'  -  ClC'C  (21) 


where 


C 


1 


ue 

2u 


(22) 
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The  momentum  thickness  8  is  a  constant  and  is  determined  from 


OD 

TtfV  -  2n  J  ux(U  -  drx 


(23a) 


0  -  ^  J  u(l  -  u)  dC 


(23b) 


It  is  related  to  the  drag  of  the  body,  D,  by 


2  2 

D  -  pU  n8 


Substituting  Eq.  (21)  for  f  '  into  Eq.  (20)  yields 


(23c) 


CfJ  +  (C  +  DfJ  +  2f2'  +  \  C^e“2C  +  CXA(C  -  l)e"C  -  0  (24) 


The  function  f 2( 0  roust  satisfy  the  boundary  conditions 


f2(0)  -  0  yC  f2(0  -  0  as  C  -  0 


f2(C)  -  0  as  C  - 


To  solve  Eq.  (24)  we  set 


f2'< C)  -  (C  -  l)e‘C  g2(C) 


*nd  g  (£)  then  satisfies  the  equation 


A  + 


2  -  C2e‘C 

; .ct--. .o.  g>  .  .  .i.... 

C(C  -  1)  S2  2C(C  -  1) 


.  - 


t 
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This  is  immediately  integrable  and  the  solution  is 


g2(C) 


-  cf  PC  “t  -C  rC  et  eC 

-T  JffVdt+  —  +c2  — 


+  CjA  c  -  J  .  ! j 


where  C_  and  a  are  arbitrary  constants.  Hence,  £2(0  is  given  by 


f2(0 


(C  -  l)e'C  V  dt  +  e'2C  +  °2  (C  •  1)e"C  £  T  dt  -  l 


+  C1A[ 


(C  -  l)e-t  in  £  -  2e 


Since 


pC  e-t  e  £ 

-  dt  ~  -  — —  +  constant 

,1.  t  t 


rc  efc  e£ 

J  dt  ~  +  constant 


the  only  term  in  Eq.  (26)  which  does  not  exhibit  exponential  decay  to 
zero  as  Q  -*  00  is  the  integral  multiplied  by  C2»  Hence,  we  set  C2  *  0. 
For  C  -  0 


e-t 

-  dt  ~  in  C  +  constant 

'  a 


(28) 


which  means  if  f'(C)  is  not  to  become  infinite  (as  in  Q  when  C  -*  0, 
we  must  choose 
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With  these  choices  of  C2  and  A,  Eq.  (26)  reduces  to 


f2'(C) 


e_t  -C 

1)  dt  -  Xn  C  +eM2 


To  check  the  second  of  the  boundary  conditions  in  Eq.  (25),  we 
first  calculate 


*5«>  -  -i 


|  e‘C  |(C  -  2)  £  dt  -  in  {  +  (i  +  l) 


e_C  +  3  -  -■ 


Then 


lim  c»fr.  1 

£-0  f2<°  '  "8 


*  lim  fC  e  C  ^  ■  ,  ,  , 

-  2  c-ho  Ja  T  dt  ‘  in  c  +  4 


and,  using  Eq.  (28),  we  see  that 


f U(C.) 

C-K)  r2U; 


constant 


Hence 


■K  ^(O  -  0 


as  C  0 


and  the  required  boundary  condition  is  satisfied. 

The  last  boundary  condition  is  f2(0)  ■  0.  This  is  satisfied  by 
choosing  the  constant  of  integration  after  integrating  f*(C).  That 
this  can  be  done  follows  from  the  fact  that  f^(C)  has  been  made  finite 
at  C  ■  0. 

Thus  the  required  solution  for  f *(C)  is  given  by  Eq.  (29).  We 
note  that  the  lower  limit  on  the  integral  in  the  solution  is  still 
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arbitrary  after  satisfying  all  the  boundary  conditions.  This  arbitrari 
ness  introduces  a  term 


ke“^(C  -  1)  (30) 

in  f^(C)  which  cannot  be  determined.  The  reason  is  that  this  function 
satisfies  the  required  boundary  conditions  for  f^CC)  both  on  the  axis 
and  on  the  edge  of  the  wake.  If  F  "(C)  represents  the  function,  then 
the  boundary  conditions  it  must  satisfy  are 

VFf'(C>  -  0  as  c  -  0 

F  '(£)  -»  0  (exponentially)  as  C  “*  05 


The  second  of  these  is  seen  to  be  satisfied  immediately  for  F'(0  - 
e”^(£  -  1).  For  the  former,  since 

F'(C)  -  e'£  -  (C  -  De'£ 


we  see  that 


F"(0)  -  2 


so  that  it  too  is  satisfied.  (Note  that  in  terms  of  r  and  x  this 
eigenfunction  can  be  written 


c‘C  (C  -  1) 


-rZ/4@ 

e 


which  can  be  compared  to  the  corresponding  eigenfunction  for  the  two 
dimensional  flat-plate  far  wake 


e-y2/4P 
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found  by  Stewartson^  and  Crane. 

The  solution  for  the  velocity  u  is  now  given  by 


u 


-  1) 


p£  e-t 

J  V  dt  “  c 

O' 


(31) 


Integrating  Eq.  (18)  we  obtain 


xx(8)  -  A  in  8  +  B 


and  substituting  into  Eq.  (10) 

x(8)  ■  0  +  AlnP+B  +  ... 

To  invert  this,  we  proceed  as  follows: 

8  «  x  •  A  In  P  -  B  +  ... 

■  x-AjPn(x-A/nP-B+...)  -B  +  ... 
«  x  •  A  In  x  ■  B  +  ... 

Since  A  ■  C^/8,  we  obtain  finally 


8 


PS  X 


--|/n 


x  -  B 


If  8(x)  is  introduced  into  Eq.  (31),  we  obtain 


u 


1  +  £i£i  +  !ie-cis_£  +  M£i 

1  +  2  x  +  16  2  +  2 

X  X 


(32) 


(33) 


(34) 
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where  D(£)  Is  some  function  of  £.  This  does  not  yet  explicitly  dis- 

2 

play  the  x  dependence  of  u,  since  £  ■  r  /43  also  involves  x.  Substi¬ 
tuting  f3(x)  into  £  we  obtain 


c 


1  C1  in  x 

x  +  8  2 

x 


+ 


B_ 

2 

x 


+ 


(35) 


from  which  it  can  also  be  shown  that 


where  £  «  r  /4x.  Then  u  can  be  written 


u 


C1  e'C  C1  -C  ,7  n  in 

l+-2~  'T6e  <C  ' 


(36) 


Along  the  axis,  at  £  =  0,  this  reduces  to 


fi  +  in  x  „  / 1_\ 
*  16  x2  lx2) 


u  (x)  =  1  + 

o 


(37) 


Equations  (36)  and  (37)  agree  with  the  corresponding  equations  given 
in  Ref.  1. 


-13- 


REFERENCES 


Berger,  S.  A.,  The  Incompressible  Laminar  Axl symmetric  Far  Wake, 

The  RAND  Corporation,  RM- 4 8 60 -ARP A ,  July  1966. 

Van  Dyke,  M. ,  Perturbation  Methods  in  Fluid  Mechanics,  Academic 
Press,  New  York,  1964. 

Stewartson,  K. ,  "On  Asymptotic  Expansions  in  the  Theory  of  Boundary 
Layers,"  J.  Math,  and  Phys.,  Vol.  36,  1957,  pp.  173-191. 

Crane,  L.  J.,  "A  Note  on  Stewartson's  Paper  'On  Asymptotic  Expan¬ 
sions  in  the  Theory  of  Boundary  Layers'  J,  Math,  and  Phys.  , 

Vol.  38,  1959,  pp.  172-174. 


